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Abstract 


The DMSP P-7 satellite is an operational Air Force aeteoro logical 
satellite which carried a magnetometer for geophysical measurements. The 
magnetometer was located within the body of the spacecraft in the presence 
of large spacecraft fields. In addition to stray magnetic fields, the data 
hare inherent position and time inaccuracies. Algorithms were developed to 
identify and remove time varying magnetic field noise from the data. 
Techniques developed for Magsat were then modified and used to attempt 
determination of the spacecraft fields, of any rotation between the 
magnetometer axes and the spacecraft axes , and of any scale c h a n ges within 
the magnetometer itself. The corrected data were then used to attempt to 
model the geomagnetic field. This was done in combination with data from 
Magsat, from the standard magnetic observatories, from aeromagnetic and 
other survey data, and from DB-2 spacecraft magnetic field data. The 
results obtained are inconsistent and contradictory. Characterisation of 
the problem is clearest when model coefficients are compared between models 
with DMSP data only and models with other data types only. In that case 
the gll term in the DMSP models is consistently lower in magnitude by 10-25 
nT from the trend inferred from models based on other data. Possible 
causes include: First, the spacecraft fields may simply be so strong or 
variable that any correction procedures are destined to fail; Second, in 
addition to fixed or slowly varying magnetic sources, there may be a 
substantial amount of soft magnetization on the spacecraft, for which we do 
not at present know how to correct; Third, the timing and position error on 
the data furnished to Goddard Space Flight Center may result in systematic 
errors; Fourth, the automated procedure for removing noise from the data 
may be inadequate and could possibly introduce error into the resulting 
data. Future DMSP missions can be upgraded in terms of geomagnetic 
measurements by upgrading the time and position information furnished with 
the data, placing the magnetometer at the end of a boom, upgrading the 
attitude determination at the magnetometer, and increasing the accuracy of 
the magnetometer. 
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I. Introduction 


The DMSP P7 spacecraft was launched on 18 Not, 1083 into a 98.74 degree 
inclination orbit, with apogee 844 ka altitude and perigee 822 lea. (Rich, 
1984) . The priaary purpose of the spacecraft was to obtain tropospheric 
aeteorological data. However , a triaxial fluxgate aagnetoaeter was included 
on the spacecraft in order to aonitor the geophysical environaent . This 
report deals with the exaaination of these aagnetoaeter data to evaluate 
their usefulness in describing the earth’s core-produced geoaagnetic field. 

The field is assuaed to be curl free and so representable by a 
potential function in the fora of a spherical haraonic series (Langel, 1987) 

n* n 

1) V = a E E (a/r)Q + l[gn“cosa^ + h n m sina^]P n B (cos0) 

n=l a=0 
n** n 

+ a E E (r/a) ° [qn®cosa^ + s n ®sina^]P n *(cos5) 
n=l a=0 

where a is the aean radius of the earth (taken to be 8371 . 2km) , r , 6 , and f 
are the standard spherical coordinates, and the Pn*(cos0) are the Schmidt 
quasi— normalized fora of associated Legendre functions of degree n and order 
a. The suaaation involving (a/r) represents sources internal to the 
satellite orbit, whereas that involving (r/a) represents external sources 
such as the ring current. 

The magnetic field is then given by 

2) B = -W 

Theoretically, 1) and 2) hold only if n* and n** go to infinity. The 
aeasured internal B contains contributions from both the Earth’s core and 
froa its crust; n* in 1) is thus chosen so that V represents fields from the 
core but not the crust. Langel and Estes (1982) concluded that the core 
field dominates for n<13 and the crustal field for n>15. The n* chosen for 
highly accurate measurements, such as from the MAGSAT satellite, usually 
equals 13, but because the DMSP observations are less accurate, n* does not 
exceed 11 in any DMSP model. 

The DMSP F7 magnetometer was mounted on the satellite body, as opposed 
to being attached to a boom, because of spacecraft engineering constraints 
(Rich, 1984). The magnetometer, a triaxial fluxgate, was aligned with the 
spacecraft X, Y and Z axes, which are defined as follows: X is vertically 

down, Y is along- track and Z is cross-track. The three sensor units were 
built by the Schonstedt Instrument Co., Reston, VA, in the 1960 s . The 
electronics unit for the magnetometer was built by the Applied Physics 
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Laboratory of Johns Hopkins University, Laurel , MD , based on the design of 
the MAGSAT fluxgate sagnetoseter . 

The sagnetoseter acquired field seasuresents at a rate of 20 sasples 
per second. Heasurenents were in the fors of counts, with one count 
equalling 12 nano-Teslas (nT) . According to Rich (1984), the instrument was 
not intended to surrey the sain geosagnetic field, so it was not calibrated 
with high accuracy on the ground, nor recalibrated in orbit. 

Because of the close proximity of the sagnetoseter to on-board 
electronic instrusentation, its data were contaainated by non-randoa 
instruaental noise, with nagnitudes of up to several thousand nT. The 
attitude of the spacecraft was seasured to an accuracy of about 0.1 degree, 
or 380 arc-seconds. While this attitude accuracy is not as good as that 
obtained with MAGSAT, in principle it is of sufficient accuracy to enable 
seaningful rector seasuresents. In the absence of other near-Earth 
satellite aagnetic field data for this tine period, and in view of the 
success of aethods used on MAGSAT to solve for spacecraft fields (see 
section III) , it was decided to investigate the possibility of processing 
the DMSP F7 data to a stage where they nay be useful for sain field 
sodeling. 

The reaainder of this report will deal with several topics related to 
the DMSP data: 1) the processing steps and problens encountered in removing 
both randos and coherent noise fron the data, 2) the correction of data 
offsets and attitude errors using the GSFC sain field software (FIT 
progran) , 3) the utilisation of processed data in sodeling the core-produced 
geosagnetic field, and 4) the evaluation of those nodels in regards to the 
usefulness of this type of data for main field studies. 


H. Transformation of On-tape Data to Magnetic Readings. 


The magnetometer data contained on the basic DMSP data tapes received 
fron the Air Force is in the form of magnetometer counts, which must be 
converted into field values in nT in order to be useful. Data is arranged 
on the tape as a header, containing time and position information for each 
minute of operation, followed by 80 magnetometer readings (1 per second) . 

All times on the data are rounded to the nearest second. All positions on 
the original Air Force data are expressed in nautical miles, rounded down to 
the nearest nautical mile. Nautical miles were converted to kilometers (one 
nautical mile equals about 1.8 km) in subsequent data processing. 

Rounding of the time and position introduces an error of unknown 
magnitude into the data. An estimate of the possible magnitude of this 
error can be made as follows. The satellite velocity is about 7.4 
ks/second. A one second time error is therefore equivalent to an along- 
track position error of 7.4 km. Combined with the possible 1.8 km roundoff 
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a rr or in position , tht tottl along- track position arror nay bo as ouch as 
0.2 ka. These errors are in addition to those due to actual error in the 
computed epheaerides . 

Lange! (1976) has estiaated the effects of orbit error on near-earth, 
polar, satellite nagnetic field experiaents. His results are suaaarised in 
Table 1: 

Table 1 


Coaponent Maxiaua gradient (nT/ka) 


Mean of gradient magnitude (nT/km) 


Vertical Along- trk Cross-trk 


Vertical Along- trk Cross-trk 


Br 

-28.0 

-13.3 

6.8 

14.9 

5.3 

2.0 

Be 

18.0 

-6.5 

23.4 

7.6 

3.0 

1.7 


8.4 

-2.0 

23.3 

2.6 

0.7 

2.1 

B 

-28.0 

-6.1 

5.7 

18.4 

2.5 

1.8 


A position error of 1.8 ka in the vertical direction will give maxiaua 
field errors of 50.4, 32.4, and 15.12 nT in By, B$, and IU, respectively. A 
position error of 1.8 ka in the cross-track direction will give maximum 
errors of 12.24, 42.12, and 42 nT, and an error of 9.2 ka along-track will 
give maximum field errors of 123, 60, and 19 nT in Bp, B$, and B^, 
respectively. It seeas probable that the resulting errors are in the 50-150 
nT range and that they may be at least partly systematic in nature. 

20 readings per second for each of the magnetometer X, Y and Z axes 
were originally recorded by DMSP, although only readings #1 and #11 were 
written on the tapes sent to the Geology and Geomagnetism branch at Goddard 
(See Appendix A for on- tape data formats) . For the Goddard aain field 
studies, only the magnetometer reading associated with the header record was 
utilized. This preserved sufficient data density (one reading per minute) 
to fully describe the main geomagnetic field. 
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The magnetometer was calibrated prior to launch at the NASA Goddard 
Space Flight Center magnetic test chamber with the following results (Rich, 
1084): 


3) Measurement 
(nT) 

(radial) fBx) 

I I 

(along- trk) I By I 
I I 

(cross-trk) IB*J 


Calibration Matrix * 
(nT/count) 

(12.1001 -0.0056 0.0103) 

I I 

1-0.0247 12.1863 -0.01011 

I I 

l 0.0060 0.0232 12.1735J 


Measurement 

+ Bias 

(counts) 

(nT) 

fCxl 

r-0.0653) 

1 1 

1 1 

* ICy 1 

1 1 

+ 159.87331 
| | 

10*1 

139.4228J 


This equation is used to compute the magnetic field in spacecraft 
coordinates in nT, given a reading in magnetometer counts. However, this 
calibration does not take into account the field from the spacecraft, which 
adds greatly to the bias vector. This vector must be determined from in- 
flight data. Later work using the FIT program (see Field Value corrections) 
accomplished this, and re-determined the bias vector as: (89nT, 8457nT, 
-1441nT) for radial, along- track and cross-track measurements, respectively. 
These values are still somewhat approximate and require small corrections 
which will be introduced and discussed in sections VIII and IX. In 
addition, for computational ease, it was decided to redefine the spacecraft 
system to be compatible with the MAGSAT coordinate system, so that the 
spacecraft X axis is defined as cross- track, Y is radially down and Z is 
along- track. The correct transformation of DMSF magnetometer counts to 
nano-teslas in spacecraft coordinates compatible with MAGSAT is thus: 


4) Measurement = Calibration Matrix * Measurement + Bias 


(nT) 

(nT/count) 

(counts) 

(nT) 

(cross-trk) 

rBx) 
1 1 

fO.0069 0.0232 12.1735 1 
| | 

rc x ) 

i i 

r-i44ii 
1 1 

(radial) 

1 l 

1 By 1 
1 1 

= 112.1001 -0.0065 0.0193 1 * 

| | 

ICyl + 
1 1 

1 89 1 

| | 

(along-trk) 

1 1 
IB*J 

1-0.0247 12.1863 -0.0101J 

IC*J 

l 8457J 


All data discussed in the remainder of this report are assumed to have been 
processed through this equation and have units of nT. The (Bx,By,B z ) 
measurement vector in equation ((4 will henceforth have the label Bspu, 
meaning the vector is in MAGSAT spacecraft coordinates and not yet processed 
through the corrections discussed in section III. 
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TIT . Field Value Corrections. 

According to Rich (1984) , the DMSP magnetometer nay be Misaligned 
relative to the spacecraft by as auch as 0.5 degree per axis, with the 
misalignment Measured to an accuracy of about 0*1 degree* AlsO| bending of 
the spacecraft body May result in further misalignment. In addition, the 
values of the three Magnetic components may be in error by a fixed bias or 
by a Multiplying factor. The PIT program has the capability to solve for 
corrections in these parameters in conjunction with the least squares aain 
field solution. The theory of this adjustment is as follows (see also 
Estes , 1983) : 

The PIT prograa computes three types of adjustments to vector 
satellite nagnetoneter data: 1) A diagonal calibration matrix containing 

*slope* parameters , which is multiplied times the measured vector to correct 
for magnetometer drift, 2 ) a bias correction vector which is subtracted from 
the measured vector to correct for constant magnitude offsets, and 3) a 
rotation matrix which is multiplied tines the measured vector to correct for 
angular offsets of the magnetometer from ideal satellite coordinates. These 
adjustment parameters are applied to the measured uncorrected data in 
spacecraft coordinates according to the equation: 

5) Bspc = TSM * TCAL * (Bspu - bias) 

where: Bspu is the uncorrected measurement vector in spacecraft 
coordinates, as given in equation 4). 

Bspc is the corrected measurement vector in spacecraft 
coordinates . 

bias is a vector of magnetometer bias corrections in addition 
to those given in equation 4) . 

TCAL is the calibration correction matrix of slope parameters. 

TSM is the rotation correction matrix. 

The elements of bias are: (BSi, BS 2 , BS 3 ) , where BSi are component 

biases derived in the FIT program, with values derived and discussed in 
sections VIII and IX. 

TCAL has elements: fl/SLi 0 O') 

I I 

I 0 1/SL2 0 I 

I I 

l 0 0 I/SL 3 J 

where SLi are slopes derived in the FIT prograa. SLi and BSi are 
applied to the satellite X axis, SL 2 and BS 2 to satellite Y axis, and SL 3 
and BS 3 to the satellite Z axis components. 
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The elements of TSUl are based on three Buler angles ( €*, Cy , and €s) 
solved in execution of the FIT program. (Note : In the FIT program, as of 
2/28/88, c z i» denoted 62, «y is denoted 61, and e B is denoted 63.) Using 
the notation TSMij, where i is the matrix row and j the matrix column, these 

are: 


8 ) 


TSMn =* cos6 7 *cos£ s 

TSU12 = cos6y*cos£x*sin6| + sin£j*sin£x 

TSM13 = -co86y*sin6x*sin6g + sin6y*cos6x 

TSM21 ■ -sine B 

TSM 22 * cos6x*cose B 

TSM23 = -sine x *cose B 

TSM31 = -siney*cose B 

TSM32 * cosey*sin«x - siney*sine B *cosex 

TSM33 ■ cos£y4*cos£x ♦ sin6y*sin6x*sine 8 


The TSM matrix in the FIT program is derived from 3 rotation matrices 
(denoted R x , Ry, Rs) is the spacecraft coordinate system. R x is a left- 
handed rotation through the angle 6x> about the spacecraft X axis. Ry is a 
left-handed rotation about the spacecraft T axis, with angle of rotation 6y. 
R s is a right-handed rotation through the angle e B , about the spacecraft Z 
axis. 


The matrices Rx, Ry and R B are thus 


(1 

1 

0 

0 1 
l 

Ry = fcOSCy 0 

1 

sinfiy 1 
1 

10 

1 

COS6x 

1 

-sine x 1 
1 

1 0 1 
I 

0 1 
1 

(0 

sine x 

1 

COS6 x J 

1 

l-siney 0 

1 

COSCyJ 


Rs 


f coses sines 

I 

I -sine B cose z 

I 

l 0 0 


0 1 
I 

0 I 

I 

1 J 


Rotations provided by R x , Ry and R B are illustrated in Figures la), b) and 
c). 

TSM is created by rotating about the X axis first, then about Z, then 
about Y, e.g. , 

8) TSM * Ry*Rs*R x 


1. Matrices TSM and TGS are named to conform to notations used in Estes, 
1983. 
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An older version of the PIT prograa, docuaented in Estes (1983) , used a 
different order of rotation, about the Z axis first, about the new X axis 
second, and finally about the new Z axis. It is now known that this Z-X-Z 
rotation order fails to adequately resolve the first and third angles when 
they are large, and so the present order of rotation was instituted, with 
successful results. 

The relation of euler angles e x , ey, e* to "roll, pitch, yaw* notation 
is dependent on the spacecraft axis designations. For exanple, since the 
MAGSAT Z-axis is pointed in its along-track direction, 6 a is roll, 

6y (radial) is yaw, and e x (cross track) is pitch. 

The PIT prograa calculates soae field quantities in earth-fixed 
cartesian coordinates . The coordinate origin is at earth 1 s center j the X 
axis points along 0° longitude; the Y axis points along 90° meridian; and 
the Z axis points along the geographic north pole. Information is therefore 
required on the relation between the corrected spacecraft measurement 
vector, Bspc, and its analog in earth-fixed coordinates Bef, for every data 
point. This information is contained in transformation matrix TGS: 

9) Bef = TGS * Bspc 

The TGS matrix itself is an approximation computed from the formula: 


10) TGS = f-n x -r x v x 1 

I I 

I ~Dy -Ty Vy I 

I I 

l-n a -r a v»J 


where: r x = cos^cosX, ry = cos^sinX, r a = sin^ 

n x * cos^qCosXq, ny = cos^nsinXn, n a = sin^ n 

v x = n y r a - r y n a , v y = n a r x - r a n x , v a = n x r y - r x n y 


X, and X n are defined as: 

f = Geocentric latitude 
X = East longitude 

= -8.74®, the inclination of the vector normal to the orbit 

X n = X ± arccos [-tan (^n) *tan (^) ] , where + is used for a 

descending orbit (N to S) and - for an ascending orbit. 
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It should be noted that this equation is not accurate for latitudes greater 
than about 75*. DMSP data between 75* and 81.26* were therefore not 
utilised in subsequent analyses. Derivation of this restriction and of the 
TGS transformation elements themselves say be found in Appendix B. 

Combining equations #5 and #0 yields: 

11) Bef = TGS * TSM * TCAL* (Bepu - bias) 

Bspet and TGS are read or computed from input data, and TSM, TCAL and 
bias are solved for in execution of the FIT program. 


IV. Noise Sources in the Magnetic Data. 


The DMSP data were examined initially by Rich (1084) , who found three 
sources of magnetic noise. The first two are high frequency sinusoidal 
signals with periods of 0.576 and 3.456 seconds. These are caused by the 
rotating X-ray scanner, designated the SSB/S instrument, which is mounted 10 
to 15 inches from the magnetometer sensors and generates a small magnetic 
field. These high frequency noise sources are of magnitude less than about 
30 nT and are not a concern in the present study. 

The third noise source found by Rich is the operation of the satellite 
torquing coils. These are turned on for durations of about 4 minutes at 
various times throughout the DMSP mission. When the coils are on, the 
magnetic field data is offset by a constant level shift of 3000 to 14,000 
nT. This type of noise is screened out by the 'gross outlier" criterion for 
reducing the data (see section VII) . 

A fourth known noise source consists of fields in the 100 - 150 nT 
range which result from turning on transmitters and tape recorders when over 
a tracking station. It is assumed that most data so affected will be 
eliminated in the various outlier tests described in section VI. 

A noise source not discussed by Rich (1084) was discovered by examining 
orbital plots of residual data, i.e. data which have had a preliminary field 
model subtracted. These residual data show strong periodic trends with 
amplitudes of up to 70 nT, after other corrections were applied. An in- 
depth discussion of this periodic noise and its removal is found in section 
VII. 
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V. Observatory and Surrey Data used in Modeling. 

D&ta froa land-based aagnetic observatories and surveys are 
instrumental in augmenting the satellite data for modeling purposes. 
Observatories have permanent locations and are useful in measuring the field 
behavior over an extended tine period. They usually measure declination, 
plus the horisontal and vertical field components, from which X, Y, Z 
magnetic components (North, Bast, and Down in a geodetic system) may be 
derived. The observatory data used in the DMSP studies were obtained from 
the World Data Center in Boulder, CO. They are annual means of daily 
measurements, so each falls on a half year (e.g., 1979.5). Since 
observatory biases (see next section) were solved for, only those stations 
with at least 3 measured components and 3 years of data were retained. 

The time distribution of observatory data used in this study is: 


Table 2 


Tine period 

1979 - 1980 

1980 - 1981 

1982 - 1983 

1983 - 1984 

1984 - 1985 

1985 - 1986 

1986 - 1987 


# of stations 

143 

151 

155 

149 

140 

91 

1 


The geographic distribution of observatories is roughly equal to that 
displayed in Langel et.al. (1988b). 

Survey data include aeromagnetic data (Project Magnet and a Caribbean 
survey) and land data. The Project Magnet data are 3— component vector 
recordings, the Caribbean data are scalar values only, and the land data may 
be declination, inclination, horisontal, north (X), east (Y) , vertical (Z) 
or scalar readings. The aeromagnetic data were averaged over a 220 km 
distance prior to the least-squares fitting procedure, in order to reduce 
data volume, attenuate crustal anomalies, and obtain an estimate of data 
quality for data weighting in the fitting procedure. The averaging method 
is described in detail in Langel et . al . (1988b) . 

The land data were checked for quality by comparing to a preliminary 
field model, and assigned an error estimate via a geographic binning 
procedure (see Langel et. al., 1988a for description of binning process). 

The data have extensive geographic coverage; the temporal coverage is 
described in Table 3. 
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Table 3. 


Ma gnet data. Caribbean Land data 


Time 

I 

Y 

Z 

B 

D 

I 

H 

Z 

Y 

z 

B 

1980 



... 


1285 

110 

761 



521 

159 

1981 

308 

306 

306 

— 

105 

13 

88 



— 

71 

— 

1982 

436 

436 

436 

— 

132 

39 

88 

— 

— 

68 

34 

1983 

276 

276 

276 

— 

105 

— 

124 

— 

— 

44 

60 

1984 

— 

— 

— 

142 

92 

— 

71 

6 

7 

71 

3 

1985 

— 

— 

— 

— 

329 

89 

151 

— 

— 

105 

151 

Total 

1018 

1018 

1018 

142 

2048 

251 

1283 

6 

7 

880 

407 


VI. Calibration Procedure 


Fluxgate magnetometers are subject to scale and zero-level drifts due 
to changes in mechanical properties and electronic components. For MA6SAT, 
these mere solved for in flight by comparison with the on-board scalar 
magnetometer data. Since no scalar magnetometer was available for DMSP, an 
attempt was made to use data from magnetic observatories as a "standard* 
against which to calibrate the DMSP data. Calibration is accomplished by 
solving for the matrix TCAL (see equation #5) in the modeling procedure. In 
this process the "standard" data constrain the main field solution 
sufficiently to permit the TCAL solution. 

Two methods of calibration were developed, both of which use 
observatory data. The use of such data in a main field model necessitates 
taking into account the presence of crustal fields, fixed measurement 
biases, and any other field not represented by the model. Langel et. al. 
(1982) accomplished this by designating the sun of such fields Bn, and the 
measured field at the observatory B. Bn is called the "bias" at a given 
observatory. B equals Bn + Bi, where Bi is the field from the spherical 
harmonic model. 

In the first method of calibration, Bn is determined for each 
observatory prior to the spherical harmonic solution, under the assumption 
that the GSFC (12/83) model, epoch 1980, is highly accurate (Langel and 
Bates, 1985). Bn for any observatory is taken to be the measured annual 
mean nearest 1980 (generally at 1980.5) minus the field value computed from 
the GSFC(12/83) model, truncated at degree 10. (For consistency, the 
truncation level is the same as that used for models derived with DMSP 
data.) If an observatory has no data for 1980 then the bias is computed at 
a later year (1981.5 or 1982.5), and that observatory is assigned a larger 
error estimate in the subsequent spherical harmonic fit. 
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Once Ba for each observatory is thus obtained, it is assumed to remain 
constant through tine. lg is then subtracted from measured obserratory 
values around the epoch of the DMSP data being calibrated. The corrected 
observatory data are then interpolated to the average epoch of the DMSP data 
and a field model is generated. The DMSP data are assigned a 200nT weight, 
compared with 5nT for the observatories. In the spherical harmonic 
solution, satellite biases, Buler angles, and the elements of TCAL are 
allowed to vary in addition to the main field coefficients, to obtain the 
calibration constants . 

In the second calibration procedure, Bg is solved for simultaneously 
with the main field coefficients. This solution requires that data from at 
least three epochs exist for each observatory (Langel et al, 1082) . Thus 
all observatory data from 1080 through the DMSP time period are used, as 
opposed to using a single epoch of extrapolated data, as in the first 
procedure. The observatories and DMSP data are given weights based on 
residuals to previous fits (generally 25 nT for observatories and 40-80 nT 
for DMSP data) . The epoch of the model is set to 1080. Both main field and 
time terms are computed, along with the DMSP satellite biases, Euler angles 
and TCAL elements. The GSPC (12/83) information matrix is used as a priori 
information, which "anchors* the constant solution at 1080. Observatory 
biases are solved for, i.e. allowed to vary, in this procedure, versus being 
held fixed in the first procedure. 

This second procedure gives calibration results very similar to the 
first procedure. However, the second procedure is more versatile, since it 
fan be used to calibrate the DMSP data over its entire time span, whereas 
the first procedure can only be used to calibrate data from a single epoch, 
since no secular variation is included in the model. Both of the methods 
rely on the GSPC (12/83) model for determining the observatory biases, with 
the implicit assumption that it is a highly accurate model and that the 
biases remain constant over the time span 1080 - 1085. 


VII. Automated Processing Procedure for DMSP Data. 


Preliminary field modeling 

A test model was generated from DMSP data at epoch 1084.04. The gi^ 
term from this model equaled -20,000.4 nT. A model derived from observatory 
data at the same epoch yielded g]0 equal to -20,883.4. The closeness of the 
two terms suggests the apparent adequacy of the DMSP data for main field 
modeling. A calibration of the 1084.04 data, using the second procedure 
described in the previous section, was also executed with the following 
results: SL1 = 0.0055, SL2 = 0.0006, SL3 = 1.0025. The nearness of these 

values to unity again suggests that the DMSP magnetometer measured the 
magnetic field accurately for that selection of data. These results 
indicated that DMSP data might be useful for main field modeling, in spite 
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of the large spacecraft fields. All of these studies were conducted using 
only & few days of data. On the basis of these results, it was decided to 
proceed with a larger quantity of data. 

The preliminary field aodel was removed from January 14-18 DMSP data to 
create residual data. Upon examination, these orbits of data showed strong 
periodicities. A spectral decomposition of the data revealed noise sources 
with periods equal to the orbit period (100 minutes) and subharmonics of 
1/2, 1/3 and 1/4 of the orbit period. Figure 2 displays a typical residual 
orbit from this tine period. The X and Y components most clearly 
demonstrate this periodic noise. Figure 3 is the associated spectrum. 

Peaks in the spectrum display these dominant noise periods quite noticeably. 
As will be seen below, one cause of the periodic noise is the need for 
adjustment of the Euler angle values in the TSM matrix in equation 5) . The 
other causes of this periodic noise is unknown. The peak-to-peak amplitude 
of the noise is about 300 nT before Euler angle correction and about 50-70 
nT after that correction. 


Automated procedure description 

A four-step procedure was followed to remove data spikes and periodic 
noise from the DlfSP data. The steps are as follows: 1) Fit a span of DMSP 

data, covering several days, with a preliminary field aodel. Subtract the 
field aodel to get residual data. Reject data points above 75° absolute 
latitude, and reject "gross outliers”, i.e., residual data with absolute 
values greater than a specified cutoff. Fit the residual data with a spline 
function and reject points which deviate nore than 2 standard deviations 
from that function. 2) Add residual data which is not rejected back to the 
preliminary field model. Then fit a new field model to this data with epoch 
equal to the average time of that data span. Solve for constant main field 
coefficients, magnetometer angle adjustments and biases. 3) Correct the 
original data with the angle and bias solutions. Use the computed field from 
step §2 to re-create the residual data, and re-do step # 1 , i.e. reject gross 
outliers and spline outliers. 4) Fit a Fourier function, which is composed 
of the 4 dominant noise periods (25 minutes, 33 minutes, 50 minutes and 100 
minutes) in a least-squares manner to the residual data. Reject outliers 
according to the Fourier fit, using the 2 a criterion as for the spline fit. 
Then subtract the Fourier function from the data. Add the result back to 
the computed field aodel from step §2, to create the final, corrected data 
set. 


Step 4) is somewhat ad hoc. Such periodic variations could arise from 
source corrections we have either overlooked or been unable to apply. For 
example, comparison of Figures 4a and 4b shows that much of the large 
periodic oscillation results from unadjusted Euler angles. It both is more 
meaningful and reliable to correct the euler angles than to remove the 
variations via the Fourier fit. For this reason the ad hoc Fourier fit 
correction is applied last. 
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To reiterate, the four step process corrects the data in several ways: 
a) It rejects individual bad points according to a field model, b) It 
rejects points above 75* which have an inaccurate attitude determination, c) 
It corrects the data with magnetometer rotations and biases solved for in a 
least-squares fit, d) It rejects points which are outliers to spline and 
Fourier function fits, and e) It removes periodic noise from the data by 
subtracting out the Fourier function. 

Figure 4 shows the same profile from Figure 2, after it has undergone 
the data cleaning process. Host of the periodic noise is gone, and the 
major spikes and outliers have been removed. 

Figure 5 is a world-wide distribution plot of 3 days of processed DMSP 
data. Note that no data exists polewards of ±75* latitude. Also, every 
orbit contains a no-data sons approximately 15 degrees long. These points 
have been rejected according to the * gross outlier" criterion, and evidently 
are positions in the orbit where the torquing coils were turned on. 

A diagnostic of the quality of a data set and the geomagnetic model 
based on that data set is the "spectrum” from that field model. Consider 
the quantity: 

13) Rn = (n+1) E [(g n “) 2 ♦ (b n ®) 2 ] 

m=0 

R n is the mean-square value over the Earth’s surface of the magnetic field 
intensity produced by harmonics of the nth degree. Models derived from the 
MA6SAT data are considered to be the nearest to noise-free of all available 
models. The presence of noise at any degree will increase the value of Rn- 
Figure 6 shows two such spectra from DMSP models compared with that for a 
MAGSAT model (Langel and Estes, 1982). The uncorrected DMSP spectrum, from 
January 14-18 data, begins to deviate from the MAGSAT spectrum at degree 9, 
whereas the spectrum from corrected data deviates much less, and not until 
degree 10. This analysis shows that the removal of noise sources results in 
an improvement of the character of the output field models. It also shows 
that the DMSP data are considerably noisier than the MAGSAT data, even when 
corrected. 


VUI. Individual Epoch DMSP Field Models 


The correction procedure was applied to 15 sub-sets of DMSP data, each 
containing several days of data. Subset epochs ranged from January, 1984 
through November, 1985. Each data set was chosen from a magnetically quiet 
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period as determined by the world-wide Kn index. Results of step 2, which 
solves for the field model and magnetometer adjustment parameters, are 
summarised in Table 4: 

Table 4. 


Date 

81° 

gl 1 

hll 

Cx 

«y 


BIAS1 

BIAS2 

BIAS3 

yrs 

nT 

nT 

nT 

deg 

deg 

deg 

nT 

nT 

nT 

84.02 

-29895 

-1927 

5522 

-.038 

-.449 

0.005 

12.0 

2.6 

-1.7 

84.06 

-29893 

-1928 

5532 

-.065 

-.446 

0.006 

0.2 

-.5 

1.4 

84.21 

-29887 

-1935 

5523 

-.114 

-.459 

-.004 

7.8 

2.7 

-8.7 

84.34 

-29872 

-1925 

5516 

-.172 

-.457 

-.013 

11.4 

-2.0 

-12.5 

84.47 

-29860 

-1922 

5514 

-.138 

-.451 

-.016 

2.4 

-7.6 

-11.2 

84.63 

-29866 

-1932 

5503 

-.063 

-.474 

-.009 

-11.9 

-12.1 

-3.8 

84.71 

-29866 

-1927 

5505 

-.092 

-.478 

-.006 

-18.4 

-9.5 

-2.4 

85.06 

-29857 

-1918 

5496 

-.050 

-.496 

-.022 

-91.2 

-56.7 

20.0 

85.34 

-29856 

-1910 

5492 

-.129 

-.471 

-.013 

-91.2 

-67.0 

4.1 

85.46 

-29838 

-1920 

5494 

-.125 

-.467 

-.017 

-93.6 

-68.3 

1.3 

85.46 

-29843 

-1916 

5491 

-.130 

-.472 

-.020 

-86.7 

-65.9 

1.3 

85.60 

-29842 

-1915 

5495 

-.098 

-.475 

-.009 

-100.4 

-74.2 

4.6 

85.75 

-29847 

-1908 

5490 

-.074 

-.520 

0.013 

-112.6 

-72.3 

7.3 

85.82 

-29843 

-1914 

5484 

-.081 

-.511 

0.021 

-113.1 

-63.7 

12.8 

85.90 

-29832 

-1905 

5489 

-.030 

-.518 

0.032 

-110.0 

-57.5 

25.7 


Figures 7a) through 7e) are derived from Table 4. They display 
solutions for g]0, gj.1, h].l, the three Euler angles, and the three biases 
for each DMSP data set throughout time. The main field coefficients 
decrease in magnitude with time as expected from earlier models, but the 
trend is not smooth. This could indicate that the data sub-sets have 
marginal geographic distribution, or that the DMSP data are not sufficiently 
stable over time. The Euler angle solutions are fairly consistent, with €y 
(yaw) varying slowly from -.44 to -.52 degrees, 6 X (pitch) averaging about 
-.1 degrees and € a (roll) averaging about aero. The bias values show a 
noticeable break between September, 1984 and January, 1985, most strongly in 
X and Y. Biases at January, 1985 depart sharply from the previous bias 
trend in all three components. This jump is evident in the biases only, and 
its cause is uncertain. One possible explanation is that on 30 October, 
1984, the solar array panel was rotated 90°. This could result in a changed 
contribution to the bias field from the solar array since both its position 
and its total current were changed. Another, though less likely, spacecraft 
change that could contribute to the bias change is that on 7 November, 1984, 
the skew momentum wheel was reset so that it drew 100 ma less current. 

The bias values in Table 4 are part of the value of the vector 
parameter bias to be used in equation 5), i.e. they are a small time 
dependent correction to be applied in addition to the large bias values of 
equation 4) . A small further correction is derived in section IX. 
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H. Field Models Using All DMSP Data. 


Final data processing steps 

After the DMSP data sub-sets were corrected according to the euler 
angles and biases in Table 4, they were prepared for field modeling using 
the entire coabined set. Before concatenation, data froa each tiae period 
fgre first exaained according to 3— hourly Kp indices and hourly DST values. 
If the Kp index was greater than 2+ or the average DST value was greater 
than lOnT or less than -30nT, then that 3 hour section of data was rejected. 
Next, data froa each period were sorted into equal area geographic bins and 
randoaly deleted in each bin until the best approxiaation to the following 
distribution was attained: 3 vector (i.e., 3 coaponents each) points per 

bin for dip— latitudes less than 30 degrees, and 9 sca l a r points per bin for 
dip-latitudes greater than 30 degrees. This distribution is based on work 
done by Lowes and Martin (1987) , which indicates that vector data are 
required at low magnetic latitudes to avoid the a perpendicular error 
effect", while scalar data at high latitudes are more effective than vector 
measurements, because they are more accurate, especially in the presence of 
field-aligned currents. 

After binning, the 15 different time periods were merged together into 
a single file and entered into a format suitable for input into the FIT 
program. The final data set consisted of 7750 X and Y measurements, 7768 Z 
data, and 16097 B data. 

All of the DMSP data were then calibrated according to the second 
procedure in section VI. The field model for this calibration was of degree 
eleven in the constant field, first derivative (secular variation or SV) 
terms, and second derivative (secular acceleration or SA) terms. The 
terminology used herein to designate such a model is (N1,N2,N3), where N1 is 
the degree of expansion of constant spherical harmonic terms, N2 the degree 
of expansion of first derivative terms and N3 the degree of second 
derivative terms. The calibration model is thus described as (11,11,11). 
Note that N1>N2>N3 . The GSFC (12/83) model was used as a priori information, 
the epoch was 1980, and the calibration solution included satellite 
biases, Euler angles, and TCAL elements. In the solution, the Euler angles 
were negligible, the bias vector equaled (4.8nT, -O.SlnT, 0.65nT), and the 
slope solution equaled: SL1 = 0.99994, SL2 = 0.99960, SL3 — 1.00119. 

These calibration results were applied to the data before creation of 
the final field models, and these are the final corrections to be applied in 
accord with equation 5). Thus, in equation 5), the value of bias is the sum 
of this bias value and those found in Table 4; the TCAL matrix is derived 
from these values of 3L1, SL2, and SL3; and the TSM matrix is derived from 
the Euler angles of Table 4 with no additional correction. 
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Use of DB-2 data in Models 

The DESP data were combined in subsequent fits with 3-coaponent and 
scalar aagnetic field data from the Dynaaics Explorer 2 (DB-2) spacecraft. 

A complete description of these data is giren in Ridgway (1088) , but is 
suaaarised here. 10,600 rector DB-2 aagnetic observations, between 10/81 
and 1/83, were aade available to the Geologjr and Geoaagnetisa branch, 

Goddard Space Flight Center. Since the priaary aission of DB -2 was to aap 
field-aligned currents at high latitudes, the data were concentrated in 
polar regions. The data were also not well distributed teaporally, being 
concentrated at epoch 1082. 

These probleas necessitated iapleaentation of a binning procedure on 
the data and deciaation within the bins, to distribute the data aore equally 
both teaporally and geographically. This procedure cut the nuaber of data 
to 10,600 vector observations, iaproved the teaporal distribution, reduced 
data concentration at the poles, and iaproved data quality by rejecting bad 
points . 

The DB-2 aagnetoaeter was booa- mounted, in order to ainimize spacecraft 
fields. Accuracy of the spacecraft attitude deteraination, however, was 
only 0 . 5 ° - 0 . 7 ° (Lange 1 et. al., 1988), resulting in bias and rotational off- 
sets in the three aagnetic components. Attitude corrections to the data 
were empirically coaputed using the FIT software (see section III) , which 
significantly reduced residuals froa calculated field models when applied. 
The DB-2 data were also calibrated to ground observatories in the same 
Banner as the DMSP data. Results from the attitude correction and the cali- 
bration are as follows: 

Euler angles : e x = 0.220°, 6 y = -0.091°, e a = 0.411° 

Biases biasi = - 10 . 8 nT, bias) = 6.9nT, bias] = 2 . 8 nT 

Calibration slopes : SLi = 0.99928, SL 2 = 0.99976, SL 3 = 1.00011 

(For biases and slopes, along-track = 1 , radial = 2, and cross-track = 3) 

Application of these corrections improved the data fit to all 
subsequent geomagnetic field models. The data residuals, relative to the 
DB -2 field model at 1982 (excluding external field solution), are: 

ff x = 44 nT 
Oy s 55 nT 
<x a — 48 nT 
ffb = 24 nT 

These results show that the scalar magnitude data (residual ffb) is fit much 
better than the vector component data. Later models , whether using DE-2 
data alone (see Ridgway, 1988) or DE -2 plus DlfSP data (see following 
section) were generally most successful when utilizing only the scalar DE-2 
data. 
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Initial Models 


Tha first aodel using all of the DMSP data, called the MODMSDE aodel, 
was of degree (11,11,11), Meaning sain field static teras, linear tine 
terns, an d 2 Q d order tine terns up to degree 11. The a M a of MODMSDE stands 
for MAGSAT a priori infomation, froa the GSFC(12/83) aodel. The a 0 a is 
for observatory annual aeans extending froa 1970.5 through 1985.5. The a DM a 
is for DMSP data, with rector data at low latitudes and scalar Magnitudes 
abore 30° dip latitude. The a S a stands for world-wide surrey data of all 
aeasureaent types, which were weighted with a correlated weight aatrix (see 
Langel et. al., 1989). The a DE a stands for DE-2 data, dirided into rector 
and scalar data just as for the DMSP set. The aodel was created in steps, in 
the order of the letters. Also solred for were obserratory biases, three 
external field coefficients of the first order, and the 4 DST Multipliers 
(one for giO, qi^, qil and s].l) . The aodel epoch was 1980, and at 1980 the 
aodel was hearily constrained by the MAGSAT data, with the other data 
priaarily determining the temporal variation terns. 


Examination of secular rariation terms 


The secular rariation results, for MO, MODM, MODMS, and MODMSDE, all of 
degree (11,11,11) at epochs 1980, 1985, and 1990, are displayed in Table 5: 


Table 5. 


Model 

Epoch 

ftl° 

ftl 1 

fil 1 

62 ° 

te 1 

B2 1 


1980 

25.7 

8.0 

-28.0 

-19.8 

5.2 

-3.3 


1985 

21.3 

9.2 

-20.4 

-17.0 

5.4 

-18.4 


1990 

17.0 

12.4 

-12.7 

-14.1 

5.6 

-33.5 


1980 

25.2 

6.2 

-20.48 

-17.2 

3.9 

-9.9 

MODM 

1985 

28.6 

8.9 

-20.43 

-18.7 

2.9 

-17.3 


1990 

32.1 

11.7 

-20.38 

-20.2 

1.9 

-24.6 


1980 

24.4 

9.5 

-23.3 

-16.8 

4.6 

-8.7 

MODMS 

1985 

29.0 

6.3 

-18.0 

-19.0 

2.0 

-18.0 


1990 

33.7 

3.1 

-12.8 

-21.2 

-0.5 

-27.4 


1980 

25.1 

15.2 

-19.5 

-15.7 

3.7 

-10.1 

MODMSDE 

1985 

28.5 

1.8 

-21.3 

-19.9 

2.8 

-18.9 


1990 

31.9 

-11.9 

-23.2 

-24.2 

2.0 

-23.8 
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The behavior of these aodels is not consistent. For the MO aodel, $1® 
decreases with tiae, in contrast to MODM, MODMS and MODMSDE, where it 
Increases iron 1980 to 1990. For gll, the MO and MODM aodels are in good 
agreeaent, the MODMS is substantially different, and the MODMSDE shows a 
totally opposite trend. For All, both MO and MODMS increase over tiae, at 
roughly the saae rate. In contrast, the MODM All reaains constant through 
tiae, and the MODMSDE decreases. 

It was now apparent that the aodels had probleas. Besides the 
inconsistent secular variation results, it was found that the calculated 
correlations between the SV and SA coefficients were highly negative, 
approaching -1. This aeaas that the solution has strong, undesirable 
interdependences which night also contribute to the discrepant secular 
variation results. The paraaeter space was selectively scaled back to 
(11,10,5), (11,8,5), and (11,8,0) aodels. Whenever SA terns were included in 
the solution, hig h negative correlations with the SV terns invariably 
resulted. When the aodels were of degree (11,8,0), nost of the high 
correlations disappeared, indicating that this is near to the maximum 
paraaeter space which can be supported by the available data. 

The MO aodel was also examined. It also showed large correlations, 
although not as high (±0.75), and more scattered in the solution space. 

These correlations persisted even when the MO solution was cut to a (11,8,0) 
aodel, indicating that available observatory data alone are also not 
sufficient to solve for a model with this sise parameter space. 

Examination was also aade of the secular variation results from the 
M0(ll,8,0) , (11,11,0) and MODM(11,10,5) models. glO, gll and All terms for 
these models are plotted in Figure 8, along with the saae coefficients 
derived from observatory first-difference data froa 1940 through 1979. If 
the new models are valid, they should be consistent with the trend of the 
first difference plots. The glO plots generally match the first-difference 
trend for the MO models. The MODM model is soaewhat off of the trend. For 
gll, both MO (11,10,5) and MO (11,8,0) fall roughly in the expected 
coefficient range. MO (11,11,0) is somewhat out of the range, and the MODM 
aodel is definitely off of the trend. For the All plot, all of the MO 
aodels and also the MODM aodel are consistent with the trend of the first- 
difference coefficients. 

Clearly the different gll results in Figure 8b are inconsistent. 
Agreement of the M0(ll,8,0) model with the first-difference trend argues, as 
above, that the observatory data distribution is inadequate. Addition of 
the DMSP data should, in principal, greatly improve the data distribution 
and permit solution for more parameters. But clearly the MODM model is in 
serious disagreement with the first-difference trend. Adding to the 
confusion is the apparent agreement of the MODM(11,10,5) model with the 
MO (11, 11) model. 
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To this point DMSP vector data hu boon utilised. Since vector dst&^is 
aore error-prone than scalar data, it was decided to examine the possibility 
that at least some of the model problems stemmed from systematic errors in 
the DMSP vector data. 


MOSPMDB (11.8.5) and (11.11.51 scalar models 

Errors in magnetometer attitude, either random or systematic, should 
not affect the measured scalar magnitude. Accordingly, in order to 
eliminate any effects of strictly attitude errors, in the next series of 
models only scalar DMSP and DE— 2 data were used, in addition to 
observatories and survey data. The data included 13,121 DMSP data points at 
an average epoch of 1985, and 5100 DE-2 points with an average epoch equal 
to 1982.3. The satellite data were added to a base model consisting of 
MAGSAT a priori, observatory data, and survey data, called MOS. The epoch 
of both the (11,8,5) and (11,11,5) models is 1980. 

In both resulting models the secular acceleration terms again have high 
correlations with the corresponding secular variation terms, ftl®, 4l* and 
Bjl coefficients for each model are displayed in Figure 9, along with 
coefficients derived from the first-differences. (For this plot, first- 
difference coefficients were also generated from 1980-1985, to add to the 
1940-1979 data used in Figure 8.) Plotted in Figure 9 are the MO (11,8,5), 
MOS (11,11,5), and the MOSDM and MOSDMDE (11,8,5) and (11,11,5) models, at 
epochs 1980 and 1985. 

Examination of Figure 9 again focuses attention on ftll. Th® ftl® 

Ell models are all in reasonable agreement. Models of ftll also show some 
sense of internal agreement, but now with a totally unexpected trend. The 
Jjl coefficient decreases sharply but steadily from 1982 to 1985, both in 
the first-difference coefficients and in the models utilizing satellite 
data. Models with observatory and survey data [MO (11, 8, 5) and 
MOS (11, 11,5)] , which have SA terms, show the same trend, though not as 
pronounced. This is highly unexpected behavior, unless some sort of 
geomagnetic jerk is in progress. 


MOSDMDE (11.11.5.31 models 

The possibility of a jerk in the field at 1982 prompted a test of a 
(11,11,5,3) model, which provides for non-linear fluctuation of the secular 
variation terms. Only scalar DMSP and DE-2 data were used, as in the 
previous models. This model was produced despite the high correlations 
inherent in solving for an overabundance of time terms. Once MOSDM 
(11,11,5,3) and MOSDMDE (11,11,5,3) models were completed, parameters were 
also cut back to obtain a MOSDMDE (11,9,5,3) model. This proved to be very 
similar to the (11,11,5,3) models. 
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Plot* of th« secular variation coefficient* fro* these aodels present 
no consistent picture. For ftl® and $2 the aodels with satellite data are 
in disagreement with the first-difference trends. For All, 82 * sod &2^ 

the trends are siailar, but not exact. The aodels for g]l show a decreasing 
tend after 1080, but the first-difference trend differs significantly fro* 
the trend of aodels with satellite data. Examination of first-difference 
plots of observatory data showed no conclusive evidence for a "jerk” in this 
period. However this is not a conclusive result since the tiae span of 
available data froa aost observatories was too short to definitively 
deteraine the trend. Our conclusion is that the data available are not 
adequate to discern the existence of a jerk. 


Trends of the spherical haraonic coefficients 

We are left with the pussle of a possibly anoaalous behavior of the &].l 
coefficient, particularly when DMSP data are coabined with observatory and 
surface data. It is not yet clear that the DMSP data are consistent with 
the surface data, or that the DMSP data are reliable for geomagnetic field 
aodeling. 

Figure 10, which plots the giO, gjl, and hil coefficients, not the 
derivatives, derived froa the IGRF (1070,1975,1980), the DE-2 calibration 
(1982) , and the individual DMSP runs from Table 4, shows more clearly the 
nature of the problem with the aodels based on the DMSP data. For the g]0 
plot (Figure 10A) , the DMSP results fall fairly well along the linear trend 
of the IGRF models. The DB-2 calibration term, at 1982, is slightly high, 
but not significantly off -trend. The gj.1 plot (Figure 10B) reveals much 
more about the DMSP data. All of the g]l coefficients derived from 
individual DMSP fits are significantly off of the trend of the IGRF models, 
having a lower magnitude than expected. The DB-2 g]l, however, is slightly 
higher than expected indicating a discrepancy between the two data sets. 

The h]l term (Figure 10C) apparently varies in a non-linear manner, with 
both the DB-2 and the DMSP-derived coefficients falling along the general 
trend of the IGRF models. 

Figure 11 allows further examination of the behavior of the computed 
gil. Figure 11A shows the individual DMSP fits, along with the DE-2 and 
DMSP calibration models. The DMSP calibration matches the individual DMSP 
fits very well. It diverges noticeably from the DB-2 calibration. Figure 
11B displays the individual DMSP fits along with the (11,9,5,3) MOSDMDE 
model, the (11,9) MOSDMDE model, and the M0S (11,9) model. This plot shows 
that the DMSP models depart noticeably from the models based only on 
observatory and survey data (MOS model) . The DMSP data require a lower gil 
term. When combining DMSP data with observatory and repeat data, as in the 
M0SDMDB(11,9,5,3) model, the gil term must change dramatically with time to 
accommodate all data types. 
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X. Conclusions 


The DMSP body-mounted magnetometer data, recorded between 1984 and 
1986, hare been examined for their utility in geomagnetic field modeling. 

The data contain aajor sources of error, as indicated by large data spikes 
of up to several thousand nT, periodic noise of several hundred nT, and 
rotational and bias errors in the aagnetoaeter of up to 0.5 degree and one 
hundred nT. 

The data spikes are corrected by a gross-outlier cutoff relative to a 
preliminary field model and a spline fit. Periodic noise is removed by 
fitting the data residuals with a 4-component Fourier series. Magnetometer 
attitude errors are corrected by solving for rotational offsets using the 
FIT software, and the biases are similarly corrected. The total magnitude 
of the data is adjusted by calibrating it aga i nst ground observatories. 

These corrections improve the data quality, as indicated by residual 
plots and by power spectra of resultant field models. However, although the 
field models show magnitudes roughly within the realm of believability , 
closer inspection reveals basic flaws. Any models using Taylor series time 
representations give high correlations between the secular variation and 
secular acceleration parameters. This is true for both vector and field 
magnitude data. More importantly, the models derived from fits to the 
individual data sets show significant differences to the trend of previous 
IGSRF models and models derived from DE-2 data. These differences are 
manifested particularly in the gil spherical harmonic coefficient. In 
hindsight, if the plots displayed in Figure 11 had been examined early in 
the analysis process, much of the confusion over the secular variation 
results could have been avoided. 

Unless further methods are found to correct these discrepancies, it is 
concluded that the data at present are not suitable for main field modeling. 
There are several possible reasons for this. First, the spacecraft fields 
may simply be so strong or variable that any correction procedures are 
destined to fail. Second, in addition to hard magnetic sources and fixed 
amplitude variable magnetic sources (such as currents) , there may be a 
substantial amount of soft magnetization on the spacecraft. This would 
result in induced fields which would vary with time and position, and for 
which we do not at present know how to correct. Third, the timing and 
position error on the data furnished to Goddard Space Flight Center may 
result in systematic errors. Fourth, the automated procedure for removing 
noise from the data, particularly the removal of Fourier components, may be 
inadequate and could possibly introduce error into the resulting data. 

Nothing can be done about the first and second items in the present 
data. If further effort were to be expended upon the DMSP F-7 data, it 
should deal with items three and four. A possible approach would be as 
follows: a) obtain new data from the AFGL, for the same quiet day periods, 
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with no round-off in tine and position information, b) Eliminate spurious 
data by plotting all of the data, in short sections (1 or 2 orbits) , and 
hand-selecting data segments which are relatively free from variable 
spacecraft fields due to equipment being turned on or off. Hopefully this 
procedure would make the Fourier correction unnecessary. 

It is not certain that these procedures would be able to salvage the 
niSP F-7 data for main field modeling purposes. However, the exercise of 
processing the data has been instructive, despite difficulties encountered, 
and good possibilities exist for use of future data. 

If the measurement accuracy can be sufficiently upgraded, future DMSP 
data could make a significant contribution to main field geomagnetism. This 
is because the present need is for long-term operational observatory-class 
measurements. Such data would not only provide operational tools such as 
field models, navigation charts and submarine magnetic detection baselines, 
but would provide a research tool desperately needed to probe the nature of 
the geodynamo, the conductivity of the mantle, and the core/mantle momentum 
exchange mechanism. Keys to upgrading the data quality are: 1) placing the 

magnetometer on the end of a boom, 2) upgrading the attitude determination 
at the magnetometer location, and 3) increasing the accuracy of the 
magnetometer . 

It is our understanding that step 3) has already been taken, at least 
in part. A boom-mounted magnetometer would significantly reduce periodic 
noise and bias errors in the data. Combined with more accurate attitude 
determination, such an addition to the DMSP magnetic measuring capabilities 
would greatly enhance the usefulness of DMSP data in deriving accurate main 
field models. 
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XI. Appendix A: DlfSP Data Format 


I. Tape parameters 

Tape density: 6250 BPI 

Logical record length: 75 bytes 

Block sise: 1875. 

Fixed block, ASCII. 


II. Data format 

The data are arranged as one epheneris record every minute, 
followed by 60 magnetometer records (one per second) . 

A. Bphemeris record 


location 

format 

description 

1- 4 

14 

Year. 

5- 8 

14 

Day of year 


9-14 

18 

Time of epheaeris, seconds U.T. 
Altitude (nautical miles) . 

15-18 

14 

19-28 

F10.2 

Geographic latitude. 

29-38 

F10.2 

Geographic longitude. 

39-48 

F10.2 

Corrected geomagnetic latitude. 

49-58 

F10.2 

Corrected geomagnetic longitude. 

59-68 

F10.2 

Corrected geomagnetic local time 

69-71 

3x 

Blank. 

72-75 

14 

Number of data records following 


ephemeris (usually = 60) . 


B. Data record 


location 

format 

description 

1- 6 

16 

Time of data record (seconds UT) . 

7-11 

5x 

Blank. 

12-29 

316 

X, Y, Z magnetometer counts for 
first of 20 samples per second. 

30-33 

4x 

Blank. 

34-51 

316 

X, Y, Z magnetometer counts for 
eleventh of 20 samples per 



second. 

52-55 

4x 

Blank. 

56-75 

1012 

Ten data quality flags. 
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TIT . Appendix B: Derivation of DMSP Transformation Matrix, T6S. 


For every DMSP magnetic reading, known position information 
includes the geographic latitude, longitude and altitude of the 
satellite, and the inclination of the satellite orbit plane. From this 
information, an operator may be derived which transforms measurements in 
the local satellite coordinate system to an earth-fixed cartesian 
coordinate system (Originally derived by D. Chinn of Science 
Applications Research Inc.). 

Figure B1 displays a DMSP satellite orbit superimposed upon 
rotating Barth-fixed cartesian coordinates. The inclination of the 
satellite orbit plane with the earth’s equatorial plane is 98.74 
degrees. A local orthogonal coordinate system may be assigned to the 
satellite, with one axis pointing along the satellite track, another 
normal to the satellite track, and a radial axis pointing outwards from 
the earth-fixed origin. The corresponding unit vectors are v, n, and r, 
respectively . 

The unit radial vector r, in cartesian coordinates, equals: 

Bl) r ■ r x i + ryj + r z k 


where : 


B2a) 

r x = cos^cosX 

B2b) 

Tj * cosjtsinX 

B2c) 

r B = sinjt 


and where f equals the latitude, X the longitude of the satellite. 

The unit normal vector n may be described in an analogous manner, 
except that replaces ^ and X n replaces X in equations Bl and B2. 
equals a constant -8.74° anywhere on the orbit, while Xn varies 
throughout the orbit. X n is solved for, given ^ and X, by utilizing the 
fact that r * a equals zero: 

B3) cos^sinXcos^nsinXn + cos^cosXcosf> n cosXn + sin^sin^n = 0. 

Dividing by cos^ and cos^n and switching tangents to the other side 
of the equation yields: 

B4) -tan^tan^n - cosXcosX Q + sinXsinX Q . 

Utilizing the addition formula for cosines yields: 

B6) -tan^tan^n = cos (X - X n ) 


or: 
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B6) X n = X + arccos( -tanftanf*n ) (- for ascending orbit, 

+ for descending orbit) . 

Once Xn is solved for, then r and n are fully determined by B2 . ▼ 

is derived froa the relation: 


B7) ▼ = n x r 


The transformation 

froa 

cartesian 

(i,j,k) to local (r,n,v) 

coordinates is 
B8) f r 1 

thus: 

fr x 

t 7 

r.l 

f i 1 

1 1 

1 


1 

1 1 

1 u 1 

= Inx 

n y 

n B l * 

1 3 1 

1 1 

1 

1 

1 1 

l v J 

'-▼x 

T y 

▼«J 

l k J 


The MAGSAT-compatible satellite coordinates (x-cross track, y- 
radially down, a- along track) correspond to -n, -r, and ▼, respectively, 
so that the rotation from cartesian to spacecraft coordinates is thus: 


B9) 


f cross'! 

f-n x 

1 

” n y 

1 down 1 = 

I I 

1 

l-r x 

l 

” r y 

lalongJ 

^x 

v y 


-n*! 

rn 

1 

i i 

-r z l * 

iji 

1 

i i 

y z J 

UtJ 


The T6S matrix is a rotation from spacecraft to cartesian 
coordinates, and so equals the inverse (= transpose for a rotation 
matrix) of the above equation: 


fil 
1 1 
IJI 

1 i 

f-n x 

1 

_r x 

▼x^ 

1 

1 1 

1 

■ 1 -n y 

1 

-ry 

rjr ! 

1 1 

(down 1 

| 1 

UtJ 

1 

l-n a 

-r z 

V*J 

LalongJ 


_____ 

TGS 





Accuracy of derivation 

The accuracy of the TGS derivation depends strongly on the accuracy 
of the Xn value in equation B6) . A small error in longitude (X) 
determination will propagate to an equal error in Xn, but an error in 
latitude (^) will propagate to a larger error in Xn* This relation is 


- 25 - 



seen by taking the derivative of X n with respect to 

Bll) dXn/d^ » -K/( cos^cos2(^) - K2 s in2(f>) ) with K = tan(-8.74) 

The relation of dXn/d^ to f is as follows: 


TABLE B1 


* dXn/d* 

81.28* • 

80* 10.4 

78* 5.2 

78* 3.3 

75* 2.8 

If l is known to within 0.2°, for instance, then the uncertainty in 
Xn would equal 0.56, 0.87, 1.0 and 2.1 at 75°, 76*, 78°, and 80*, 
respectively. Because of this, it was decided not to accept data above 
75° latitude. 

Alternative formulation 


An alternative approach to defining the TGS aatrix could be taken 
if certain DMSP orbital parameters were available. The following 
derivation is included for clarity of the meaning of TGS. Consider a 
typical satellite orbit in fixed inertial space, e.g., Figure B2. The 
X-Y-Z coordinate system is fixed in inertial space with Z along the 
Barth’s rotation axis and X along the vernal equinox. The vernal 
equinox is the ascending node of the Earth’s orbit around the sun. The 
epoch of X and Z must be specified. Y is chosen to make a right handed 
system. For our purposes, consider the plane of the orbit to be fixed 
in the X-Y-Z system. Define the spacecraft coordinate system, X’-Y’-Z’ 
as follows: X’ will be positive in the radial direction, Z’ will be 
normal to the orbit plane, with Y’ chosen to make a right handed system. 
For circular orbits, Y’ is along track. The angle of the orbital plane 
with respect to the equatorial plane is called the inclination, i; the 
intersection of the orbital plane and the equatorial plane, at the point 
where the satellite is going north is the ascending node. The angle 
between X and the ascending node, Q, is called the right ascension of 
the ascending node. The satellite position in the orbit is then 
specified by its angle in the orbit plane, w, measured from the 
ascending node. 

To transform a vector or position from the X’ system to the X 
system requires a set of three rotations: (1) a rotation, Ri, around the 
Z' axis through the angle w; (2) a rotation, R2, through the angle i 
around the new X axis (the axis through the ascending node) : and, (3) a 
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rotation, R 3 , around the Z axis through the angle Q. This gives 



fcosw 

-sinw 

O'! 

B12) 

Rl 3 Isinw 

cosw 

01 

l 0 

0 

1J 



f 1 

0 

0 1 

B13) 

R 2 = 1 0 

cosi 

-sini 1 

l 0 

sini 

cosiJ 



fcosQ 

-sinfl 

O'! 

B14) 

R 3 = Isinfl 

cosQ 

01 

l 0 

0 

1 J 


If R = R 3 R 2 R 1 > then 

Rll 3 cosQcosw - sinflcosisinw 

Rl2 = -cosQsin* - sinOcosicosw 

Rl3 = sinflsini 

R21 = sinQcosw + cosQcosisinw 

R22 = -sinQsinw + cosQcosicosw 

R23 = -cosQsini 

R31 = sinisinw 

R32 = coswsini 

R 33 = cosi 


If the X-Y-Z coordinate system is the true-of-date geocentric 
inertial system, called Celestial Coordinates, CC, then the 
transformation to Barth fixed coordinates is as follows. The rotation 
of the Earth, L, since the reference time of the coordinate system is 


BIS) L = DFRADY*6. 3003881 + DELTDY*0 . 1720279 + GHA 

where DELTDY is the number of days since the reference time of 
the coordinate system, 


DFRADY is the fraction of day 


GHA 


is the Greenwich hour angle at the reference time 
of the coordinate system. 
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If the X-Y-Z coordinate rector is denoted X and the Barth-fixed 
coordinates, X«, Y« , Z«, are denoted X», then X« is giren by 


B16) U « TLX 


where 


r cos(L 4 X) sin(L ♦ X) O') 

B17) Tl * I -sin(L + X) cos(L + X) 01 

10 0 1 J 


where X is the longitude. 

To then transfer froa the Earth fixed equatorial system, X«, to the 
topocentric N-B-V (North, East, Vertical, spherical Barth, geocentric) 
systee, the transformation matrix is giren by 


B18) 


Tt - 


f cos (l + r/2) 


I 0 

l-sin(^ + r/2) 


0 sin + r/2) ") 

1 0 I 
0 cos (l + r/2) J 


where 9 is the geocentric colatitude and ^ = r/2 - 9 is the latitude. 


Combining results, we then have 

f crossl 

B19) X« = TLBX’ = QX’ = Q’l down I 

l alongj 


where Q = TlR and, in the notation of equations BO) and BIO) , X* is the 
rector (-down, along, -cross) , and it is easy to show that 

Q’ll = -Qi 3 = -cos(L + X)Ei 3 - sin(L + X)R23 

Q’l2 = -Oil = -cos(L + X)Bn - sin(L + X)R21 

Q*13 * Ql2 = cos(L + X)Ri2 + sin(L + X)R22 

0*21 = -023 = sin(L + X)Ri3 - cos(L + X)R23 

Q’22 = -021 = sin(L + X)Rn - cos(L + X)R21 

Q*23 = 022 = -sin(L + X)Ri2 + cos(L + X)R22 

0*31 = -033 = -R33 

O’ 32 = -031 = -R31 

0’33 = 032 = R32 
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Comparison of B19) with BIO) shows that TGS = Q* or that 
rx = Qli; Ty = 021; r, * Q31; n x - 0l3; n y = 023; “s = 033 
and that 

B20) tanX n * n y /n x = 023/013 = cos(Q+L+X)/sin(Q-L-X) . 
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TABLES 


Table # Description Page 

1 Effects of orbit error on near-earth, polar, 3 

satellite aagnetic data, for radial, theta, 

and phi coaponents and scalar field magnitude. 

2 Temporal distribution of observatory annual 9 

means data used in DMSP models. 

3 Temporal distribution of land and aeroaagnetic 10 

data used in DUSP models. 

4 glO, gjl, and hil coefficients, plus Euler 14 

angle and bias solutions, generated from fits 

to 15 sets of DMSP data spanning 1984 to 1986. 

Table 4 data is plotted in Figure 7. 

5 ftl?, ftl 1 , hi 1 * $1°, El 1 17 

Si 1 terms for MO, MODM, MODMS, and 

M0DMSDE(11 , 11 , 11) models. 

B1 The relative error of X n > u approximated by 26 

equation B6, as a function of latitude of the 
DMSP satellite. 
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Figure Captions 


Figure # 


Caption 


1 Relation of Euler angles to spacecraft coordinate system, 
describing rotations about the X, Y, and Z axes, respectively. 
Rotations around X and Y (figures a and b) are left-handed, 
whereas Z is right-handed. 

2 DMSP orbital X, Y, and Z magnetic component data which have 
had an estimated field model removed, revealing strong 
periodicities in the residuals. The dashed line is a spline 
fit to the residuals. 

3 Power spectra of X and Y DMSP residual data from Figure 2. 

4 Y-component of DMSP data from Figure 2, demonstrating removal 
of outliers, magnetometer rotation and bias correction, and 
subtraction of Fourier periodic function. The dashed line is a 
spline fit to the residuals. 

5 World-wide geographic distribution of typical three day 
sequence of processed DMSP data. Note the absence of data 
above 75° latitude and the large mid-latitude data gaps where 
the torquing coils were activated. 

6 Geomagnetic field spectra. Rq is the total mean square 
contribution to the vector field from all harmonics of degree 
n. The spectra are from field models produced from 
uncorrected and corrected DMSP data; the Magsat spectrum is 
included for comparison. 


7 Plots of g]0, gjl, hil, Euler angles, and biases versus time 
(yrs) , for field model solutions from 15 DMSP data sets 
spanning 1984 - 1986. 

8 Plots of fti°, ftl 1 , and hi* for M0(ll,8,0), MO (11, 11,0), 
M0(ll,10,5) and MODM(11,10,5) models, at epochs 1980 and 
1985. Coefficients derived from observatory first 
differences from 1940 through 1979 are included for 
comparison . 
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0 


Plots of &0, ftll, and 111* for M0(ll,8,5), M0S(11,11,5) , and 
both MOSDM and U0SDUDE(11,8,5) and (11,11,5) models, at 
epochs 1980 and 1985. First difference coefficients from 
1940 through 1985 are included for comparison. 

10 Plots of giO, gil and hxl coefficients for individual DUSP 
fits, the IGRF 1970, 1975 and 1980 models, and the DE-2 1982 
calibration model. 

11 Plots of the gil term for individual DIISP fits. Figure 11A 
includes the DUSP and DE-2 calibration models for comparison. 
Figure 11B includes the U0DUSDE(11,9) and (11,9,5,3) models, 
and the U0S(11,9) model for comparison. 

B1 Relation of DUSP orbit with local a, r, v coordinates to 
earth-fixed cartesian coordinates. 

B2 Geometry of satellite orbit in inertial space. X-Y-Z is 

coordinate system with Z along the Earth’s rotation axis, X 
along the vernal equinox, and Y so as to give a right handed 
system. X’-Y’-Z’ is the spacecraft coordinate system with X’ 
radial, Z’ normal to the orbit, and Y’ so as to give a right 
handed system. Y’ is along track. The inclination is i, Q 
is the right ascention of the ascending node, and u is the 
angular position of the spacecraft measured from the 
ascending node. 
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Figure 3 
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DISTRIBUTION OF DMSP DATA: 
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